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Abstract 



It has been suggested that the origin of cosmic rays above the GZK limit 
might be explained by the decay of particles, X, with mass of the order of 
10^^ GeV. Generation of heavy particles from inflationary quantum fluctuations 
is a prime candidate for the origin of the decaying X particles. It has also been 
suggested that the problem of non-singular galactic halos might be explained 
if dark matter originates non-thermally from the decay of particles, Y, such 
that there is a free-streaming length of the order of O.lMpc. Here we explore 
the possibility that quantum fluctuations might account for the Y particles as 
well as the X particles. For the case of non-thermal WIMP dark matter with 
unsuppressed weak interactions we find that there is a general problem with 
deuterium photo-dissociation, disfavouring WIMP dark matter candidates. For 
the case of more general dark matter particles, which may have little or no 
interaction with conventional matter, we discuss the conditions under which X 
and Y scalars or fermions can account for non-thermal dark matter and cosmic 
rays. For the case where X and Y scalars are simultaneously produced, we show 
that galactic halos are likely to have a dynamically significant component of X 



scalar cold dark matter in addition to the dominant non-thermal dark matter 
component. 
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1 Introduction 



There have been a number of surprising observations made recently in cosmology. One 
has been the discrepency between simulations of galaxy halo formation and small- 
scale structure based on cold dark matter (CDM) and observations, which indicate 
that there is more small-scale structure than observed and that the halos are much 
less singular than predicted [^. In particular, numerical simulations of CDM halo 
formation show that the abundance of dark matter subhalos and so dwarf galaxies 
within a galaxy should be the same as the abundance of subhalos within galactic 
clusters, whereas observations indicate that the abundance within galaxies is much 
less [0]. Various explainations have been put forward, including self-interacting dark 
matter ^, halos made of condensates of ultra-light scalar particles and warm 
dark matter 0. Recently it has also been suggested that dark matter could originate 
non-thermally from heavy particle (or topological defect) decay 0. 

Another puzzle is the existence of ultra-high energy cosmic rays (UHECR) beyond 
the 10^^ GeV GZK cut-off [|], |^. A simple possible explaination of the UHECR obser- 
vations would be to have a heavy decaying particle of mass ~ 10^^ GeV and lifetime 
longer than the age of the Universe, such that high energy proton, nuclei and photon 
primaries can originate within the 50Mpc mean free path for energies greater than the 
GZK cut-off A natural origin for such heavy decaying particles is from quantum 
fluctuations generated during inflation Jll] . 



In this paper we consider whether the decaying particles which may explain non- 
singular halos and those which account for UHECR could both originate from quantum 
fluctuations during inflation. The halo problem is then assumed to be solved by a 
density of non-thermal dark matter (NTDM) coming from the decay of heavy particles 
which we label Y, such that the non-singular halos are a result of a free-streaming 
length ^ O.lMpc 0, whilst the UHECR are explained by decaying particles X with 
masses of the order of 10^^ GeV 0. 

The paper is organized as follows. In section 2 we discuss the general conditions 
for solving the galactic halo problem via NTDM from heavy particle decay. We also 
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briefly review the conditions required to account for UHECR via heavy particle decay. 
In section 3 we consider the special case of weakly interacting massive particle (WIMP) 
dark matter with unsuppressed weak interactions, in particular the conditions required 
to evade energy loss of the WIMPs via scattering from the thermal background whilst 
preserving primordial nucleosynthesis. In Section 4 we discuss whether the galactic 
halo problem and UHECR can be explained by X and Y particles generated during 
inflation and the possible observable consequences of this. In Section 5 we present our 
conclusions. In an Appendix we discuss the scattering cross-section and energy loss 
per scattering for a Majorana fermion scattering from thermal background particles. 

2 General Conditions for NTDM and UHECR from 
Decaying Particles 

2.1 Non-Thermal Dark Matter 

The problem of singular galactic halos and excessive structure on small-scales can be 
solved if the dark matter has a free-streaming length of the order of O.lMpc 0]. To 
achieve this in the case of heavy dark matter particles, the dark matter must have a 
non-thermal origin, such as a decaying massive particle^. Suppose the dark matter 
particles originate from decay of a massive particle at t = tii. Then the co-moving 
free-streaming length is given by 



A..= P^<ia«^^^(l+ta(^)) , (1) 

Jaa a n a^qUEQ \ \cinr/ / 

where is the scale factor at matter-radiation equality, is the scale factor at td, 
V is the dark matter particle velocity, and we have used the fact that v oc a^^ once the 
freely propagating dark matter particles are non-relativistic. If we assume that the 
initial dark matter particle energy is /?my, we obtain 

(^EgHEQ V \ rc JJ 



^An alternative is to have warm dark matter, a thermal distribution of particles of mass approxi- 
mately 1 KeV §. 



where 

To = (3) 

and is the dark matter particle mass. With a = 1 at present, is equal to 
the present velocity of the dark matter particles. With geq = 4.3 x 10~^(f2m/i^)~^, 
Af5 = O.lMpc is obtained for M 



X 



r, ^ 2.4 X 10-*^ ( "^^ 1 . (4) 

This is the condition for dark matter particles from Y decay to account for non-singular 

J?') this then fixes the Y mass as a function of 

the decay temperature, 

rnv « 4 X 10^ ri^U^l) (,^) GeV . (5) 



10-^/3 J \giT^) J VlOO GeV 

where g{T) is the number of degrees of freedom in thermal equilibrium and is the 
photon temperature at present. 

In addition, the dark matter particle density, assumed to be dominated by x, must 
satisfy ~ 0.3. (We assume throughout a flat Universe with a cosmological constant 
Q\ ~ 0.7.) Assuming that the number of dark matter particles from each Y decay is 
e, we have 

Pc 

where pc is the critical energy density and ?T,y(Ty) = g(T^)T^ / g{Td)T^ is the Y number 
density at present in the absence of Y decays. The requirement VL^ ^ 0.3 will impose 
a constraint on the reheating temperature after inflation, Tr. 

2.2 UHECR 

The observed UHECR are not correlated with conventional sources of cosmic rays and 
are consistent with an isotropic distribution |p. This is consistent with the idea that 
they originate from decaying massive particles [0. The particles must have a mass 
mx ~ 10^^ GeV and lifetime tx ~ 10^^(^x/3 x 10""^) yr, where ^x ~ ^x/^x 
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fraction of halo CDM in the form of X particles. For masses 10^^ GeV ^ uix ^ 10^^ GeV 
there is no conflict with the diffuse 7-ray background or positron flux in cosmic rays 
[P]. Therefore, as far as X particles generated during inflation are concerned, the only 
restriction is that VLx < 0.3. 



3 Non-Thermal WIMP Dark Matter 

The conditions in Section 2 are quite general and apply to dark matter particles 
regardless of their interaction with particles in the thermal background. For the more 
specific case of Y particles decaying to WIMP dark matter and to a significant number 
of charged particles or photons, we must impose two other conditions. 

The first condition is that the WIMPs from Y decay do not lose too much energy 
by scattering off of thermal background particles, since it is assumed that the WIMPs 
evolve by freely expanding after Y decay. 

The condition for the WIMPs not to scatter is that, for the case of relativistic 
WIMPs, the scattering rate for WIMPs from thermal background particles should sat- 
isfy = n(T)asc^E < EH, where n{T) = ^''^^^^'^ is the number density of particles 



in the thermal background E is the energy of the WIMP and AE is the energy 
lost by the WIMP per scattering with a particle in the thermal background. In the 
following we will consider WIMPs with masses of the order of myy and unsuppressed 
weak interactions, such that the scattering with thermal background electrons and 
neutrinos is via t-channel Z° exchange. Experimentally, Dirac fermions and scalars 
are excluded as dark matter particle in this case [|T2|, so we consider the case of Ma- 
jorana fermion WIMPs. In the Appendix we give the scattering cross-section for the 
scattering of relativistic Majorana fermions from thermal background particles and 
discuss the energy transfer per scattering. 

The scattering cross-section times energy lost per scattering integrated over centre- 
of-mass scattering angle for Majorana WIMPs is (Appendix) 

ascAE ^ — K = [log [—^j - 4 1 ; ^^xEr >m\z. (7) 
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and 

32na'^k'EtE^ ,40 
^scAE^ y \ K = ^ ; 4E^E^<ml,. (8) 

Here a = g^gT-fAir ^ 0.006 (with and gr the couphngs of x ^-^id the thermal 
background particles to Z" (Appendix)) and E^ ^ 3T is the energy of the thermal 
background particles. 

For the case where AE^E^ > m\ ^ (Td > Max(m|, m2)/(3/?my)), no-scattering 
imposes the constraint 

SuPkT, \PJ\0.006J \100J \10j ' 

where fc^ = ( 45^"^^ ) ^ • This is ruled out if F particles are to be generated by quantum 
fluctuations during inflation, since particles with my > Hj ^ 10^^ GeV cannot be 
generated during inflation |ll| , where Hj is the expansion rate during inflation f\. 

For iE^Er < m\^z i'^d < Min(m|, m^)/(3/3my)), we obtain an upper bound on 
the decay temperature 

T < ( ^ V^' k'Jfmzm^ 

Numerically this gives 

5.2x10-3/0.006^/2/ \( mz \ /lOO GeV\=^/V 10 V'^ 

-pyl llOO GeV j V90 Gevj V my ) \g(Tdj J 

(11) 

(We keep mz explicit in this in order to check the effect of increasing the mass scale 
of the gauge interaction.) 

Combined with the free-streaming constraint Eq. (^, we find that in general no- 
scattering requires that 

"^--'°-n^J [-p)\W-^) liooGivj IgoGivj [-^) 

(12) 



^It may be possible to generate such heavy particles via preheating after inflation ml 
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and 

(13) 

Thus in order to both evade scattering and act as a source for non-thermal WIMP 
dark matter, my ^ lOTeV and ^ 2 x 10~^GeV is necessary. 

In the above we have imposed the condition that the WIMPs effectively lose 
no energy through scattering with thermal background particles. However, if < 
Min(m|, m^)/(3/5my) then the rate of loss of energy decreases as the particles lose 
energy, since Tgc^E/E oc E^^. So in this case, if is larger than the upper bound in 
Eq. (|TB|), it may appear that it is possible for the WIMPs to lose energy and stop scat- 
tering from the thermal background whilst still remaining non-thermal, so allowing 
to evade the upper bound Eq. ([13|) from no-scattering. However, it is easy to see that 
this is not the case: the effect of WIMPs losing energy to the thermal background is 
to effectively reduce the initial energy /5my at until the no-scattering condition is 
satisfied. However, the upper bound from no-scattering and free-streaming, Eq. ([I3|), 
is P independent, so reducing P will not allow to evade the upper bound. 

A second condition, which applies if ^ T„„c ~ IMeV and if there is a significant 
number of photons produced in the Y decay cascade, is that the decaying Y particles 
should not dissociate light elements produced during nucleosynthesis |]T5[. We will 



impose the constraint following from the non-dissociation of deuterium D, using the 
conservative bound based on assuming that a fraction of the energy from Y decays 
is entirely in the form of threshold 2.3MeV photons, which is the maximum number 
possible from the cascade decay ||T6|. D non-dissociation then requires that 



n2.3p(2.3 MeV) < , (14) 

where 71,2.3 = f-yPriTd) /2.3MeV is the number density of 2.3MeV photons, n^, is the 
number density of D nuclei and p{E) is the probability of a photon of energy E 
photo-dissociating a D nucleus [l^ 

P{E) = ^ . (15) 
6 



an is the D photo-dissociation cross-section, ax is the total photon scattering cross- 
section from background particles (for threshold photons these are given hj an = 3mb 
and (Tt = 125mb 0]) and Ue is the electron density. Using Ue = ub = riBs(T), 
where is the baryon number density, rjB is the baryon number to entropy ratio and 
s{T) = ^iP"^ entropy density, this results in a T^-independent bound on the 

Y mass 

my < 2.3 — MeV . (16) 

45 (Td f-y^lxPc 

Numerically this gives 

<noa ^ f \ [ 0.3\ / 



where critical density is pc = 7.5 x lO""^^/;,^ GeV^. However, since em^ < uiy, this 
results in the condition 

/, < (^^) . (18) 



~ /i2 V5 X lo-iiy V^x/ 

Therefore unless only a very small fraction of the Y mass ends up in electromagnetic 
final states, deuterium dissociation rules out Y decay after nucleosynthesis. 

Taken in conjunction with the free-streaming and no-scattering constraint, which 
requires that T^^ ^ 2 x 10^'^ GeV, we see that in most cases the deuterium constraint 
rules out WIMPs with unsuppressed weak interactions as non-thermal dark matter. 
If we wish to retain WIMPs with unsuppressed weak interactions as dark matter 
we must have Y particles decaying dominantly to WIMPs with little or no electro- 
magnetic states produced in the decay. This is difficult, since in general a neutral 
weakly interacting particle will come in an SU{2)l representation together with elec- 
tromagnetically charged particles. So we would require a model in which the mass 
of the Majorana WIMP x" lighter than its charged SU{2)i partners such 
that 2m^o < my < 2m^±. In this case Y — > 2x° would be kinematically allowed but 
Y — i> x~^X~ disallowed. Since the mass splitting of the SU{2)l representation must 
come from SU{2)l breaking and so must be at most of the order of mw, the Y mass 
must also be of the order of mw in any NTDM WIMP model with unsuppressed weak 
interactions. 
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Alternatively we can consider dark matter particles which interact more weakly 
with the thermal background than WIMPs with unsuppressed weak interactions and 
raise the upper bound above Tnuc- This can be done either by reducing the strength 
of the WIMP coupling to the exchanged particle or by increasing the mass scale of 
the exchanged particle. From Eq. (p!3| ) we see that this requires that (m^/a^/^) is 
increased by a factor of about 25 to have an upper bound on greater than IMeV, 
where we consider rriz to represent the mass scale of a general exchanged particle. 
From the free- streaming condition Eq. we see that if ^ IMeV then the decaying 
Y particle must have a mass greater than 40TeV for of the order of 100 GeV, or 
more generally that niY/m^ ^ 400. 

In []TB| it has been suggested that dark matter neutralinos could act as non-thermal 
dark matter. They find that for the case of neutralinos with unsuppressed Z° couplings 
(Wino limit), NTDM is ruled out, in agreement with the results derived here, but that 
for neutralinos with suppressed Z° couplings (Bino limit), which scatter from the 
thermal background via slepton exchange, and with heavy sleptons of mass of the 
order of ITeV, it is marginally possible to have non-thermal neutralino dark matter 
with a free-streaming length of O.lMpc. 

4 Particle Densities From Quantum Fluctuations 

In the previous sections we have considered the conditions for the dark matter particles 
from Y decay to account for non-thermal dark matter. We now consider the origin of 
the heavy particles which may be responsible for NTDM and UHECR. In the following 
we will consider x Y particles without imposing restrictions coming from their 
interactions with conventional matter. For example, it is possible that the non-thermal 
dark matter particles and decaying Y particles could belong to a hidden sector which 
couples only very weakly if at all to conventional matter. In this case could take 
any value without disrupting nucleosynthesis. 

We will make the simplest assumptions regarding inflation, namely that the expan- 
sion rate H is fixed during inflation, with inflation followed by a coherently oscillating 
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inflaton matter dominated period characterized by a reheating temperature Tr. 



4.1 X and Y Scalars 

We will consider the amplitude of the quantum fluctuations to be such that only the 
mass term in the scalar potential plays a role i.e. we neglect self-interactions. We also 
consider the case where = at the beginning of inflation (where represents the 
X oi Y scalars), so that there is no significant contribution to the number density of 
scalars from an initial classical expectation value for cf). For m < Hj, the scalars may 
be considered massless during inflation. 

During each interval 6t ^ Hy^, the scalar field receives a quantum fluctuation 
on horizon scales 6(j) ~ Hi/2tt, which is thereafter stretched beyond the horizon as 
a classical fluctuation. Therefore after N e-foldings of inflation, the mean squared 
magnitude of the classical scalar field as seen on sub-horizon scales will be |Tl[ 

..2 NH] 



An"' 



(19) 



as a result of the random walk of the scalar field due to quantum fluctuations. Once 
inflation ends and the Universe enters the inflaton matter dominated regime, the super- 
horizon wavelength modes begin to re-enter the horizon. We denote the expansion rate 
when the fluctuation of wavelength A re-enters the horizon as if a- The evolution of 
the modes then depends on whether they re-enter before or after H = m. The modes 
obey the equation of motion 

+ 3HS^ - k^6(f) = -m^S(j) , (20) 

where k is the wavenumber of the mode. Modes entering at Hx < m have S(j) constant 
until H = m, after which they oscillate about 50 = with 5(f) oc Xjc?!'^. Modes 
entering at Il\ > m are additionally suppressed, since their equation of motion is 
initially dominated by the term. As a result their energy density evolves as 50 oc 1/a 
until k^ oc a^^ < m^, after which the mass term dominates and the modes evolve as 
50 oc The average number density of particles at if < m will therefore be 
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given by 

n ^ m < 5(i)'^ >= n^H^<m) + n(^H^>m) = {Nt - N^)m5(l)l + m < Scpfjj^^^^ > , (21) 

where < ... > denotes spatial average, 50(/f^>m) is the total contribution of modes 
with Hx > m, is the number of e-foldings before the end of inflation at which 
modes re-entering at if = m leave the horizon, Nj- is the total number of e-foldings 
of inflation (A^^ ^ 65) and 6(j)c is the amplitude of each Hx < m mode. In general, 
a mode re-entering at H during inflaton matter domination exits the horizon during 
inflation at 

N{H) = hn (^) . (22) 

With Hi ^ 10^3 GeV and i/ = m > 100 GeV, for example, A^^^ < 8. Thus to a 
reasonable approximation we can ignore the contribution of modes with Hx > m and 
simply consider A^^, = in Eq. (|2T|). In this case the gradient terms in the 6(f) equation 
of motion can be neglected and all modes simply oscillate about the minimum of the 
potential such that 5(f) oc a^'^/^ once H < m, whilst remaining constant aX H > m. 
Therefore during matter domination by the inflaton {H oc a~'^/^), (50c is given by 

50e ^ -f^ . (23) 

Thus during inflaton matter domination, the number density of (f) scalars is 

<H) « . (24) 

The number density at temperatures T < Tr is then 
4.2 X and Y Fermions 

Quantum production of fermions during inflation is possible only if the fermions are 
massive [jl^. The resulting density during matter domination has been estimated to 
be 

n = C^mH'^ , (26) 
10 



where Cq, ~ 3 x 10^'^ for a transition from inflation to matter domination. In general, 
the density of scalars will be much larger than that of fermions of the same mass, by 
a factor 

' ^ ^ > 1 . (27) 



As a result, if one has scalars and fermions with similar decay rates to dark matter 
particles (for example, supersymmetric partners), then the number of dark matter 
particles produced will be determined by the scalars. 

4.3 NTDM from Scalars 

In order to account for NTDM from the decay of Y scalars we require that Vt^ = 
en{T.y)'m^/ Pc = 0.3. Using Eq. (|25|) we find that this requires that 

^ ^ em^NT7THjg{T,)T^ ' ^ 

Numerically this gives 

/60\//i2\/ my \ AOO GeV\ /I0i3 GeV\ V^xA ^ / x 

4.4 UHECR from Scalars 

The requirement that Qx < 0.3 implies the mx-independent constraint 

Tn < 0.3 rr2 ■ 30 



Numerically this gives 



r,;280.= f™)V|^)GeV. (31) 



Thus a low reheating temperature is necessary in order to allow UHECR to be ex- 
plained by scalars generated by quantum fluctuations. 
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4.5 Relationship Between X and x Densities from Scalars 

It is important to note that the mass density of X and Y scalars from quantum 
fluctuations is approximately the same, since the number density Eq. (pS]) is inversely 
proportional to m. (This is in contrast with the case of fermions, where the number 
density is proportional to m.) Therefore, in general, the present density of dark matter 
particles from Y scalar decay is related to the density of X scalars by 

^^h'^, (32) 

where /q, as discussed below, parameterizes the uncertainty in the number density of 
scalars from quantum fluctuations. Since em^ < my, in order to have Vt^ > Qx and 
so have dark matter primarily in the form of non-thermal x dark matter, we must 
have fn > 1. If > Qx then my/e > > my / fo.^, and so for e and not very 
different from 1 and my must be of the same order of magnitude. 

parameterizes the fact that the number density Eq. (^) has an uncertainty. 
For example, if the total number of e-foldings of inflation Nt is much larger than 65, 
then most of the number density in Y scalars would be due to superhorizon modes, 
resulting in a coherently oscillating field 5(f){t) in the Universe. The initial amplitude 
of this field is determined by a 1-dimensional random walk of the field during inflation 
due to quantum fluctuations. For n steps of size ±a, the probability of a particle being 
at X is described by the normal distribution 

e 2a^n dX 

The root mean square value is given by xl^^ = no?. In our case x = 50, a = Hj /{2n) 
and n = Nt- As an estimate of the range of values /q can reasonably take we 
integrate Eq. (^) and exclude the range of values (0, and {\6(f>2\,oc) for which 

the probability is less than 0.316, chosen to give the range of values of fu which have 
a probability of 0.8. This gives = O.41|50rms| and |502| = |50rms|- (Note that the 
best estimate of the mean value of \6(f)\ is somewhat smaller than |50rms|- If we take the 
value for which the probability of being larger or smaller is 0.5, then the mean value is 
0.7\6(f)rm,s\-) The number density n oc can therefore reasonably take values in the 
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range 0.34 to 2.04 times the mean density. So in this case fn = riy/nx can reasonably 
be in the range 0.17 to 5.9, with the probabihty of being smaller or larger than these 
values being 0.316^ = 0.1. (This is independent of the definition of the mean density.) 
So long as /n > 1, the dark matter density can be primarily due to non-thermal x 
dark matter. However, we still expect to find a significant contribution to the total 
dark matter density from X scalars. This is important, since in this case galactic halos 
will be composed of a combination of free streaming non-thermal dark matter particles 
and a smaller but possibly dynamically significant component of conventional cold dark 
matter due to X scalars. This may result in different predictions for halo formation 
and small-scale structure as compared with the limiting cases of pure non-thermal 
dark matter or pure cold dark matter. 

4.6 NTDM from Fermions 

In this case, in order to account for = 0.3 we require that 

AiT^g (T^) CoT'^em^mY 

Numerically this gives, 

. 4.3 X 10- ^ (15^) (^) GeV . (35) 

Since my ^ Hj ^ 10^^ GeV for fermions generated from quantum fluctuations, there 
is a lower bound on the reheating temperature as a function of my. 

For example, we consider two mass scales for x and Y particles which are of partic- 
ular interest: the weak scale mw and the mass scale 10^^^^^ GeV associated with X 
particles and Hj. Since the largest possible reheating temperature after inflation is 
Tr ^ {HiMpiY/^ ^ 10^*^ GeV, the case ~ mw is only marginally compatible with 
Y fermions from quantum fluctuations and only if my ^ 10^^ GeV. The free-streaming 
constraint, Eq. (1), then implies that the Y fermions decay at temperature 

T. . 2.5 X 10^ C-^] (m f {,^\ GeV . (37) 



13 



(In this case WIMPs are ruled out as dark matter by the no-scattering constraint, 
Eq. (P).) For the case of very large x mass, ~ 10^'^ GeV, the reheating temperature 
must satisfy 

.....x.O=^(i^)(-^)(^?|)oeV, ,3B, 

where we have imposed my 10^^ GeV since my > em^, whilst the Y fermions decay 
at 

(In this case WIMPs are ruled out as dark matter by the nucleosynthesis constraint.) 

4.7 UHECR from Fermions 

In this case the constraint Vtx < 0.3 implies that 

Numerically this gives 

Thus a wide range of reheating temperatures is consistent with UHECR from X 
fermions. 



5 Conclusions 

We have considered the conditions under which non-singular galactic halos and UHECR 
can be explained by decaying particles produced during inflation by quantum fluctu- 
ations. For the case of WIMP non-thermal dark matter from decaying particles, the 
requirement that WIMPs with unsuppressed weak interactions do not lose energy by 
scattering from the thermal background combined with the requirement that their 
free-streaming length is of the order of O.lMpc implies that the Y particles decay at 
temperatures less than that of nucleosynthesis. Thus unless Y particles can decay to 

14 



WIMPs without producing a significant number of pliotons in tlie cascade, WIMPs 
witli unsuppressed weak interactions are ruled out by pliotodissociation of deuterium. 
It may be marginally possible to suppress the weak interactions of the WIMPs suf- 
ficiently that the Y particles can decay above 1 MeV, for example in the case of 
supersymmetry with Bino dark matter and with heavy sleptons of mass ~ ITeV [|T^. 

An alternative possibility is to consider the case where the dark matter particles and 
Y particles belong to a sector interacting only very weakly if at all with conventional 
matter, such that there is no danger to nucleosynthesis from Y decay products. In 
the case where UHECR and NTDM originate simultaneously from X and Y scalars 
generated by quantum fluctuations, the present mass density of X scalars will naturally 
be of the same order of magnitude as that of the non-thermal dark matter. In this case 
we expect the dark matter in the galactic halo to be a combination of free-streaming 
non-thermal dark matter and conventional X scalar cold dark matter, which could 
alter the predictions for galactic halo and small-scale structure formation from the 
case of pure non-thermal dark matter. This would allow an observational test of the 
idea that decaying scalar particles generated during inflation can explain both UHECR 
and non-singular galactic halos. 



Appendix. Majorana Fermion Scattering Cross-Section 



The scattering cross-section in the CM frame is 

/ da\ 1 \pf\ 



\M\' , (42) 



\dVt) ^j^j 647r2s I Pi 

where / are the initial and final three-momenta and s is the CM energy squared. 
We consider a Majorana fermion x of mass scattering from a massless thermal 
background fermion r. For simplicity we will consider a head-on collision of x with 
r along the x-direction; this can generally be arranged by boosting the frame in the 
cases where the collision is at an angle, resulting in shifts in the energy and momenta 
by factors of the order of 1. The weak interactions are described by 

Cint = ^Xl'l5XZ, + grTYil - l^)rZ, . (43) 
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In numerical estimates we use for and g-r the value of the vjjjj^Z" coupling from 
the Standard Model, g^ = gx/'^ = 92/^Cos9w = 0.19, where g2 is the SU{2)l gauge 
coupling. The spin-averaged amplitude squared for t-channel Z° exchange is then 

IGo^o'^ 

\M\^ = j^rZ^^[{Pi-<li){P2.q2) + (Pi-g2)(P2.gi) +mj(gi.g2)] , (44) 

where pi,2 arc the initial and final x 4-momenta, qi^2 ai'c the initial and final r 4- 
momenta and = {q2 — qiY- For ^ E^-, we find 

1fifl2 '2 CZ74 4 2 2 

- (^2! ^2)2 ^2 [l+2r^(l+Co.g)+2r;^(l+Co.^g) + ^^(l-Co.g)] , (45) 

where 6 is the scattering angle in the CM frame, 77 and 7 are defined by 

£;2 ^2 ^2 
r? = ^2 ^ ^ .4gN 



and 



-8£;2£;2(i-Cos^) 



Also s = + AE^Et-. The energy loss per scattering in the LAB frame (the rest 
frame of the thermal background) for a Majorana fermion scattering at angle 9 in the 
CM frame from a thermal background particle, in the limit E^ ^ m^, E^-, is given by 

2ElEr{l-Cose) , , 

^E = — ^—^ . f48) 

In the hmit AE^Er > ^ the cross-section times energy lost per scattering 
integrated over scattering angle is then 

In the limit AE^Er < m^^ ^ ^^'^ cross-section times energy lost per scattering is 

a.eA£; ^^Sf, ; = 40/3. (50) 
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